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Abstract. Classical orthogonal polynomials in one variable can be charac- 
terized as the only orthogonal polynomials satisfying a Rodrigues formula. In 
this paper, using the second kind Kronecker power of a matrix, a Rodrigues 
formula is introduced for classical orthogonal polynomials in two variables. 



1. Introduction 

One of the most important characterizations for classical orthogonal polynomials 
in one variable (Hermite, Laguerre, Jacobi and Bessel) is the so-called Rodrigues 
formula (see, for instance, 0)- 

Using this kind of formula we can write the n-th classical orthogonal polynomial 
in terms of a n—th order derivative. In fact, if we denote by {P„}„ a classical family 
of orthogonal polynomials in one variable, then 

k rl n 

(1) Pn{x) = ^" {<p{x) n UJ ( x )), n = 0,1,2,..., 

where k n is a constant, 4>{x) is a polynomial of degree less than or equal to 2, 
independent of n, and u)(x) is an integrable function in a appropriate support set. 

If degcf) — 0, Hermite polynomials appear, up to a linear change in the variable. 
If degcj) = 1, Laguerre polynomials are obtained, and if deg</> = 2, we can deduce 
two families of polynomials, Jacobi polynomials when <f)(x) has two simple roots, 
and Bessel polynomials when <p(x) has a double root. 

Formula Q is called Rodrigues formula, honoring B. O. Rodrigues who estab- 
lished the formula in 1814 for Legendre polynomials. 

Orthogonal polynomials in two variables which are solutions of partial differential 
equations were systematically studied by H. L. Krall and I. M. Sheffer (JI])> m 1967. 
They defined classical orthogonal polynomials in two variables as the sequences of 
orthogonal polynomials {Ph,k}h,k>o such that every polynomial Ph.k, with h + k = 
n, satisfies the second order PDE 

(2) L[w] ee aw xx + 2bw xy + cw yy + dw x + ew y — X n w, 

where a(x,y) = ax 2 + d\x + e\y + fx] b(x,y) = axy + d 2 x + e 2 y + fa; c(x,y) = 
ay 2 + d 3 x + e 3 y + / 3 ; d(x,y) = gx + hi, e(x,y) = gy + h 2 , and A„ = an(n- \)+gn. 
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The special shape of the polynomials involved in the above equation is a direct 
consequence of the fact that every orthogonal polynomial of total degree n must 
satisfy the same PDE. Krall and Sheffer showed that, up to a linear change in the 
variables, there are nine different sets of orthogonal polynomials satisfying such 
type of PDE. 

The first reference to a Rodrigues formula for classical orthogonal polynomials in 
two variables appears in the classical text by P. Appell and J. Kampe de Feriet (PQ). 
Later, P. K. Suetin (JHDj an d Y. J. Kim, K. H. Kwon and J. K. Lee (^01) consider 
an analogue of the Rodrigues formula for Krall and Sheffer classical orthogonal 
polynomials in two variables. In fact, for n a positive integer, they define 

(3) Pn-iA*> v) = — d *~ l d l(p n ~ l <i 

where w(x, y) is a weight function over a simply connected domain, and a symmetry 
factor of L, the linear differential operator defined in l|2[l. and p(x,y), q(x,y) are 
polynomials related with the polynomial coefficients in (J2J. Then, under some ad- 
ditional hypothesis, © defines an algebraic polynomial in two variables orthogonal 
to all polynomials of lower degree (see example 3, in Section 

The above Rodrigues formula runs only for classical orthogonal polynomials as- 
sociated with a positive definite moment functional, since it needs a weight function. 
Nevertheless, H. L. Krall and I. M. Sheffer founded classical orthogonal polynomials 
in two variables associated with a non positive definite moment functional which 
has a symmetry factor (see L. L. Littlejohn JS]), but not a Rodrigues formula like 

© (HDD- 

On the other hand, tensor product of two classical orthogonal polynomials in 
one variable, defined by 

Ph,k{x,y) = Rh(x)Sk(y), h,k>0, 

where {Rh}h>o and {Sk}k>o are Hermite, Laguerre, Jacobi or Bessel polynomials, 
satisfies a Rodrigues formula as In fact, Ph.k{x, y) can be written as a product 
of the respective Rodrigues formulas 

Ph,k{*> V) = — d * w i) — d v (*a h, k > 0. 

However, tensor products of classical orthogonal polynomials in one variable are 
not classical according to the Krall and Sheffer definition since they do not satisfy 
equation J2J, except for Hermite and Laguerre polynomials. 

Recently, the authors (see extended the concept of classical orthogonal 

polynomials in two variables to a wider framework, which, of course, includes the 
Krall and Sheffer definition and tensor products of classical orthogonal polynomials 
in one variable. 

The vector representation for orthogonal polynomials introduced in |12l 113] . and 
developed in |17| is the key to introduce the concept of classical orthogonal poly- 
nomials in two variables. Let {P n }„ denote a weak orthogonal polynomial sequence 
(see Section 2), it will be called classical (in an extended sense) if there exist non 
singular matrices A„ £ _A4 n +i(R), such that, 



(4) 



£pP* ] = div ($VP;,) + f *VP* = P* t A. 
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where 

$=( a b \ d-a x -b y \ 
\ b c J ' * \e-b x -cy )> 

and a, 6, c are polynomials in two variables of total degree less than or equal to 2, and 
d, e are polynomials in two variables of total degree less than or equal to 1 , and div 
and V denote the usual divergence and gradient operators in two variables. Observe 
that the left hand side of (@J generalizes the left hand side of the Krall and Sheffcr 
PDE 10, without any restrictions on the polynomial coefficients. Moreover, this 
new definition also include the tensor product of classical orthogonal polynomials 
in one variable. In the Krall and Sheffer case the matrices A n are scalar matrices, 
and in the tensor product case, they are diagonal non-singular matrices. 

In this paper, we will obtain a matrix Rodrigues type formula for classical orthog- 
onal polynomials in extended sense. Denoting by <J>{ n } the second kind Kronecker 
power of the matrix $ (see Bellman 0), we will show, under some hypothesis, that 
the expression, 

(5) Q* = — div {n} ($ {n} u;), n > 0, 

to 

provides a classical WOPS, where oj(x, y) is a symmetry factor of the PDE and 
div {n} is a n-th order differential operator. 

This formula generalizes in a natural way the Rodrigues formula proved in |16j , 
and JUj , and the Rodrigues formula for tensor product of classical orthogonal poly- 
nomials in one variable. 

Moreover, using our results, we will deduce a matrix Rodrigues formula for 
classical orthogonal polynomials associated with a non positive definite moment 
functional whose PDE has a symmetry factor (see example 6, in Section H3) . 

The structure of the paper is as follows. In Section 2 we collect the necessary 
basic tools. Section 3 and 4, are devoted to introduce classical orthogonal polyno- 
mials in two variables and symmetry factors associated with the partial differential 
equation The matrix Rodrigues formula JSJ, as well as some examples are 
studied in Section 6, and finally, the proof of the main result is given in the last 
section. 

2. Orthogonal polynomials in two variables 

First, we introduce some notations. Let V denote the linear space of real poly- 
nomials in two variables, and V n the subspace of polynomials of total degree not 
greater than n. 

Let A4h X k0$.) and Mhxki'P) denote the linear spaces of hx k real and polynomial 
matrices, respectively. When h — k, the second index will be omitted. 

Let A be a matrix, we denote by A 1 its transpose, and by det(A) its determinant. 
As usual, we say that A is non-singular if det(A) ^ 0. Furthermore, we introduce 
Ih as the identity matrix of dimension h. 

Moreover, we define the degree of a matrix of polynomials A 6 Mhy.k{V), as 

deg^4 = max{degOi,j(x, y), 1 < i < h, 1 < j < k} > 0, 
where di t j(x,y) denotes the (i, j)-entry of A. 
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Before discussing our approach, we briefly give some general properties and tools 
about bivariate orthogonal polynomials. For an exhaustive description of this and 
another related subjects see, for instance, [1191 HUIITT1 IT21 O UTA HT] . 

Let {nh,k}h,k>o be a double indexed sequence of real numbers, and let u : V — > K 
be a functional defined by means of the moments Hh,k — ( u i % y k ), h, k = 0, 1, 2, . . ., 
and extended by linearity. Then, we will say that u is a moment junctional. 

Some elementary properties about moment functionals acting over polynomial 
matrices A e Mhxk and B e Mhxi are given by (see pfl M HH1 ITTj ). 

(1) (it, A) = ((u,a hj ))^ =1 G M hxk {^), where A = (ai,j)^f =1 , 

(2) (Au,B) = (u,A* B). 

We say that a polynomial p(x, y) S V n is orthogonal with respect to u if 

{u : pq) = 0, VqeV, degq<degp. 

Then, we can define 

V n = {pe V n /(u,pq) = Q,Vq £ V n -i}. 

A moment functional u is called quasi definite if dim V n = n + 1. 

Definition 2.1 (H|). A polynomial system (PS) is a vector sequence {FVi}n>o such 
that 

Pn = (-P«,0j ■ • • j -fo,ra)* G -^ / (( n+1 ) xl (P n ), 

where {P n ,o, Pn— • • • > -Po,n} are polynomials of total degree n independent mod- 
ulus V n -\. 

Observe that a PS is a sequence of vectors whose dimension and total degree 
are increasing: Po is a constant, Pi is a column vector of dimension 2 of bivariate 
polynomials of total degree 1, P2 is a column vector of dimension 3 whose elements 
are bivariate polynomials of total degree 2, and so on. 

Definition 2.2 ([HI)- We will say that a PS {P n }n>o is a weak orthogonal polyno- 
mial system (WOPS) with respect to a moment functional u if 

( M ,P„P^> - 0, n^m, 
(u,P„P*J = H n , n = 0,l,2,... 

where H n 6 A / ( Tl +i(K) is a non-singular matrix. 

In the particular case where H n is a diagonal matrix, we will say that the WOPS 
{P„}„>o is an orthogonal polynomial system (OPS). Moreover, if H n = I n +\, we 
call {P n }n>o an orthonormal polynomial system. A moment functional u is quasi 
definite if and only if there exists a WOPS with respect to u (jHl)- 

In addition, a WOPS is called a monic WOPS if every polynomial contains only 
one monic term of higher degree, that is, 

P Kk {x,y) =x h y k + R(x,y), h + k = n, 

where R(x,y) G V n -\- And finally, we have that for a quasi definite moment 
functional u, there exists a unique monic WOPS associated with u. 

In this paper, we will need some differentiation tools. In fact, we will use the gra- 
dient operator V, and the divergence operator div, defined as usual. The extension 
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of this operators for matrices is introduced in 0E10IH1- Let A, -Bo, -Bi G Mhxk{V) 
be polynomial matrices. We define 



VA 



d x A 
d v A 



e M 2 hxk(P), div 



B 



d x B + d y Bi£Mhxk{V). 



We extend these definitions for n > 1. In fact, writing = V, and div^ = div, 
if we denote £>™ = (™)<9™ _4 <9*, i = 0,1,..., n, we can introduce the differential 

operators and div^ by means of 



(6) 
(7) 



V*"U= (V$A,V? A,--- ,VIA) 1 e M {(n+1)h)xk (V), 

n 

diy {n} (B ,B lr -- ,B n y = Y,V? Bi e M hxk (V), 

i=0 



where A, Bq, B\, . . . , B n G M.hxki'P) are polynomial matrices. In addition, we 
establish V™A = A, and div {0} A = A. 

The previous definitions can be translated to the linear space of moment func- 
tional using duality. For n > 0, we define the n-th distributional gradient operator 
and the n-th distributional divergence operator acting over moment functionals in 
the following way 

fp \ [po\ 

Pi i x Pi 

} = (-!)'>, div {n} 



(div 



{»} 



\PnJ 

( U °\ 

Ui 



\Un) 



> = (-l)«5>,£>f ft ), 



\PnJ 



,P) = (-!)"< 



( u \ 

Ml 



i=0 



n 



M n} p) = (-ir52(ui,v? P ). 

i=0 



3. Classical orthogonal polynomials in two variables 

In El HI Ej , the authors extended the concept of classical bivariate orthogonal 
polynomials. In fact, we define classical orthogonal polynomial starting from a 
matrix partial differential equation with matrix coefficients, a direct generalization 
of the partial differential equation studied by H. L. Krall and I. M. Sheffer in |14j . 
and P. K. Suetin in [THj . 

Let a{x,y), b(x,y), c{x,y) be polynomials in two variables of total degree less 
than or equal to 2, and let d(x, y), e(x, y) be polynomials in two variables of total 
degree less than or equal to 1. Define the partial differential operator L acting over 
V by means of 

L[p] = ad xx p + 2bd xy p + cd yy p + dd x p + ed y p. 

The operator L[] preserves the degree of the polynomials, that is, L[P n ] C V n . 

Observe that L is the left hand side of the partial differential equation (0) studied 
in |14| and |16| . without restrictions on the polynomial coefficients. 

Define 

d^ 



(8) 



$ = 



a b 
b c 



e M 2 (V 2 ), 



*= r e^xi(Pi). 
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Then, we can write 

L\p] = div ($Vp) + &Vp, 

where * = * - (div $)*. 

Definition 3.1. Let u be a quasi definite moment functional, and let {P n }n>o be 
the monic WOPS with respect to u. Then, we say that u is classical (or {P„},i>o is 
a classical WOPS), if there exist non-singular matrices A„ G Ai n+ i(M.), such that, 
for n > 1, 

(9) L[P<J ee div ($VP^ + fc'W^ = P* n A„, 
and 

(10) det(u,$)^0. 

Differential equation can be write also for n = 0, taking Ao = 0. 
Observe that, for n=l, then VP* = I2, and equation © can be written as 

div $ + * f = P^Ai =^ **=P^A X . 

The non-singular character of Ai implies that d and e are independent polynomials 
of exact degree 1. 

In ^2] an d US] 1 every orthogonal polynomial of total degree n satisfies the same 
partial differential equation, and therefore A„ = A„ I n +i, is a scalar matrix. 
Let L* be the formal Lagrange adjoint of L, defined by means of 

(11) L*[u] ee div ($V«) - div (#w), 

then it satisfies (L*[u],p) = (u,L\p\), Vp G V. 

Observe that, if u is classical, then L*[u] = 0. In fact, for any n > 0, 

F*) = (u,LK}) = ( U ,P*„) A n = 0, 

since Ao = 0, and (it, P*J = 0, for n > 1, using the orthogonality of the polynomials. 

In the authors obtained several characterizations for bivariate classical or- 
thogonal polynomials. In particular, it is shown that, under some mild regularity 
conditions, a quasi definite moment functional u is classical if and only if u satisfies 
the matrix Pearson-equation 

(12) div ($ u) = u, 

where $ and W are the same polynomial matrices defined in (J^J. Equivalently, u is 
classical if and only if 

(13) $Vu = *u. 

4. The symmetry factor 

As L. L. Littlejohn did in |15j . in order to obtain a Rodrigues formula, we consider 
symmetry factors for the differential operator £[•]. 

Definition 4.1 (QSl). We say that L[-] is symmetric if L[-] = L*[-]. L[-] is sym- 
metrizable if there exists a nontrivial function ui(x,y) such that it is C 2 in some 
open set, and is symmetric. In this case, u> is called a symmetry factor for 

L[-]. 
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As it is well known, if u) is a symmetry factor for L[-], then the operator ujL can 
be written in a compact form, 

(14) uL[.\ =div [w$V.]. 

From this equation, we can show that a function w is a symmetry factor for L[-} 
if and only if uo satisfies the matrix Pearson-type equation 11311 . 

Proposition 4.1. Let to be a nontrivial function. Then, lo is a symmetry factor 
for equation 0), if and only if u) satisfies the matrix Pearson-type equation 

(15) $Vcj = *w. 
Proof. From equation JHJ, we get 

uj LpP*] = w[div($ VP^J + $'V PJJ = 

= div[w $ V P*J - (V V P* t + w** V P*, = 
= div[w$VP^] - ($Va;-*w)*VP^, 
and the result follows using (|14|) . □ 

Observe that (|15fl is equivalent to 

(16) div ($ w) = w. 

The explicit expression for equations (|15fl and (|16|) are 

(aw) 2 + (few) y = dw, 1 f auj x + bujy = (d — a x — b y )u), 

(bu>) x + (cu>)y = eui, J \ frcj^ + ct^ = (e - 6 X - c y ) u>. 

Obviously, the nontrivial solutions of the above system of partial differential 
equations give us the symmetry factors for L[-]. 

Proposition 4.2 ( 10 ). Suposse that ac — b 2 ^ 0. Then, the differential operator 
L[-] is symmetrizable if and only if 



0_ 



c(d -Ox — by) - b(e - b x - c y ) 



b 2 



a(e -b x - c y ) - b(d - a x - b y 



b 2 



d_ 

d x 

In ^U], the authors shows that in the Krall and Sheffer case, the existence of a 
symmetry factor is a necessary condition for the existence of an OPS solution of 
the PDE 10 • A quite similar proof shows that the result is also true in the general 
case. 

Proposition 4.3. If the differential equation 0j has an OPS as solutions, and 
ac — b 2 ^ 0, then L[-] must be symmetrizable. 

5. The second kind Kronecker power of a matrix 

The second kind Kronecker power is defined in |2J p. 236, for a square 
matrix A = (ciij) of dimension 2. Let z — At be the linear transformation defined 
by A, that is, 



(18) 



z\ \ _ / &o,o a o,i \ ( t\ \ _^ | z i — a o,at\ + a>o,ifa, 
Z2 ) \ ai,o oi.i I \ fa ) \ zi = ainh + aii fa. 



Each of the (n + 1) homogeneous products z™~ l z\, i — 0,1, ... ,n, is transformed 
under l|18|l into a linear combination of the (n + 1) homogeneous products t™~ l t % 2 , 
i = 0, 1, ... Then, the square matrix of order (n + 1) specifying these linear 
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transformation is known as the second kind n th Kronecker power of A, which 
we denote by A^. 

Throughout, we will adopt the convenion A{°> = 1. Observe that = A. 

Moreover, we can check that 



l 0,0 



2 ao,o do; 



dO, 0-1,0 ( a o,o a i,i + a o,i a i,o) a 0,l 
2 CL\fi fli, 1 



«1,0 



ah 



It is possible to give the explicit expression for the entries of A^ n \ n > 0. In fact, 
if we denote A^ — (aj™^)™j =0 , then, a direct calculation shows that 



(19) 



>} 



£ 

k=0 



n — i 
k 



I 

j - k 



a, 



n — i — k k i—j+k n j — k 



0,0 



°0,1 a l,0 



1,1 ' 



< i, j < n, 



where, as usual, (™) = 0, if m < I. 

From the above explicit expression we can obtain recurrence formulas for the 
second kind Kronecker power of a matrix A. 

Lemma 5.1. There exist two recurrence formulas for A^ n \ n > 1. In fact, for 
< j < n, we get 
Recurrence I: 



Recurrence II: 



{n} 

%/ 
M 
l n,j 



Z 0,i 
An} 



a , a 



{n-l} 

hi 



+ a ,i a t 



{n-l} 



< i < n- 1, 



(n-l} , {n-l} 
fl l,0Vl,j + «1,1 a n-lj-l- 



{n-l} . {"-1} 
a 0,0 a o,j + a 0,l a 0,j-l ' 

{n — 1} . {n— 1} -, ~ ■ . 

| . +ai,ia>_ 1)J -i 1 , 1 < z < n. 



Ira a matrix form, we can express the above formulas as follows 
Recurrence I: 



A {n} 



( ao,o 

o 



o \ 

Qq.o 



V o 

Recurrence II: 
/ Qq.o 

A {n} = 



( a 0,l 

o 



ai,o / 







«0,1 



V o 

/ «o,i 



A {n-l}jl 

A ^n-l- 



ai,i ) 







ai,o 



V o 



o 



ai,o / 



A^Ll_ x + 



ai,i 



V o 



o 

Ol,l / 



where L^ l _ 1 , k = 0, 1, are n x (n + 1) matrices defined as 



(20) 



/ 1 



V o 



o \ 

1 it) 



and L n _ x 



/ 1 



V o 



a^-Vl^, 



\ 

1 / 
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6. The matrix Rodrigues-type formula 

In this Section, we shall assume that the differential operator L[-\ is symmetriz- 
able, and we denote by a; a nontrivial symmetry factor of L[-\. Our main result is 
stated in the following theorem. 

Theorem 6.1. Let u be a classical moment functional, and let 0) be the matrix 
partial differential equation associated withu. Letu be a nontrivial symmetry factor 
of L[-\, and let us assume that there exist polynomial matrices ^/oi^i G M-i^P\), 
such that 

(a$) x + (&$)„ = $* , 

(21) (b*) x + (c$) v = *9i. 
Then, for n > 0, the expression 

(22) Q*, = — div {n} ($^> w), n > 0, 

to 

provides a 1 x (n + 1) polynomial vector of degree n, such that 

(u,Q n p) = 0, VpeV n -i. 

Moreover, if {Q„}„ is a PS, then it is a WOPS associated with u. 

As we will show later, we can establish sufficient conditions in order to obtain 
WOPS from the matrix Rodrigues formula. 

Corollary 6.2. In the hypothesis of Th,eorem, \6.1\ if $ is a diagonal matrix, then 
{Qn}n * s a WOPS associated with u. 

Corollary 6.3. In the hypothesis of Theorem \6.1\ if the {n + 1) square matrix 

(uM n} ), 

is non-singular, n > 0, then {Q n }„ is a WOPS associated with u. 

Next, we will study some particular cases. Moreover, we will deduce that the 
Rodrigues-type formula given by P. K. Suetin f|16|) and revisited in Y. J. Kim et 
al. (513) is a particular case of the matrix Rodrigues type formula (|22|l . 

Example 1: The diagonal case 

Let us assume that $ is diagonal, i. e., b = 0. In this case, the matrix Pearson- 
type equation l(T7|l reduces to 

(23) (aw) I = rfw ) (ccj) y = eu>, 

and the 2x2 matrices VP?; = (ipij)j t j = o, k = 0, 1, in l(2"T|) are polynomials of degree 
less than or equal to 1 satisfying 

(a$) x = $* , (c$) y = $*i. 

This condition holds if we take V'o l = no = 0, "0o o = 2 a x , f/'i i = 2 c y , and lt 
4>o,oi such that 

a x c + a c x — cipi i, a y c + ac y = aip Q0 . 
In this way, if, for instance, the expressions 
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are polynomials of total degree less than or equal to 1, then (|21|) holds. This is 
possible, for example, when a and c are equal up to a multiplicative constant factor, 
or a y = c x — 0, that is, the polynomial a only depends on x, and c only depends 
on y. This last situation corresponds to the case where the moment functional is a 
tensor product of univariate moment functionals. 

Moreover, the second kind n-th order Kronecker power of is again a diagonal 
matrix whose elements are given by 

<^ } =a"-V, z = 0,l,...,n, 

and the Rodrigues formula gives 

Qn = (Qn,0, Qn-1,1, ■ ■ ■ j Qo,nY , 

where 

In addition, from Corollary 16. 21 we deduce that {Q„}„ is a WOPS associated to u. 

Example 2: Tensor product of classical orthogonal polynomials in one 
variable 

Tensor product of two families of classical orthogonal polynomials in one variable 
(Hermite, Laguerre, Jacobi or Bessel), {Rh}h>o and {<Sfc}fc>o is defined by means 
of 

Ph,k( x >y) = Rh(x)S k (y), h,k>0. 
These families of two-dimensional polynomials are classical according to our defini- 
tion (see IB]). The symmetry factor for the partial differential equation is given 

by 

uj(x,y) = uj 1 (x)uj 2 (y), 
where uji(x) and u>2{y) are symmetry factors of the differential equations for the 
polynomials {Rh(x)}h>o and {Sk(y)}k>o, respectively. This is a particular case of 
Example 1, since 

and 

(awi) x = du>i, (cw 2 ) 9 = ew 2 , 
are the respective Pearson equations for the univariate polynomials. Observe that 
a y = d y = 0, that is, a and d only depend on x; and c x = e x = 0, that is, c and e 
only depend on y. Then, 121(1 holds and the Rodrigues formula gives 

= {Qn,0, Qn-1,1, ■ ■ ■ , Qo,nY , 

where 

Qn- i<i (x,y)= f^—d^ia^^—dKc 1 ^), 0<i<n. 

Observe that this expression provides Q n ^i_i(x, y) as the product of the one dimen- 
sional Rodrigues formulas for {Rh}h>o and {Sk}k>o (see, for example, 0), up to 
a constant factor. 

General tensor product of univariate classical orthogonal polynomials is not in- 
cluded in the Krall and Sheffer ([H]), and Suetin (JSj) classifications for the bi- 
variate classical orthogonal polynomials, except for Hermite and Laguerre families. 
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Example 3: The Suetin Rodrigues formula 

Under several conditions, P. K. Suetin (^Sl) and Y. J. Kim et al (JOj) proved 
a Rodrigues formula for some classical orthogonal polynomials in two variables in 
the Krall and Sheffer sense. We will show that this Rodrigues formula is again a 
particular case of the matrix Rodrigues formula given in (122t . 

By elementary algebraic manipulations, equation Ijl7|l can be reduces to the form 

fa 0\foj x \_f(3\ f a uj x = (3 u, 

\ a J \ uj y J \ 7 J \ auj y = >yu, 

where a = det (3 — c(d — a x — b y ) — b(e — b x — c y ), and 7 = — b(d — a x — b y ) + 
a(e — b x — c y ). 

In the Suetin case, the special elections for the polynomials a, 6, c, d and e, 
makes deg(a) < 3, deg(/3) < 2, and deg(7) < 2. Then, if we assume 

a = ai<22, b — aibiCi, c = c\C2, 

where a± 7^ 0, c\ 7^ 0. Then 

a = a 1 c 1 a , = Ci/3 , 7 = 0170, 

and to satisfies the following system of partial differential equations 

poj x = f3 u;, gw 9 = 7 w, 

where p — a\ ao, and q — c\ ao. 

In an d the authors prove that if [a,\) y — [c\) x — 0, deg(p), deg(g) < 2, 
and deg(/?o), deg(7o) < 1, then the sequence {P n } n , with 

Pfl = (Pn,0, Pn-1,1, • ■ ■ 7 Po,nY , 

defined by means of 

to " 

provides a WOPS satisfying the partial differential equation @. 

This is again a particular case of Example 1, since Suetin hypothesis (ai) y = 
{c\) x = implies l|2"T|) . Matrix Rodrigues formula ({2"2"|l provides, for n > 0, a WOPS 

Qn = (Qn,0: Qn-1,1) • • • > Qo,n) j 

where 

Q„_ vi (x, y) = -L Q fl£-< ^ (p»-< ?* w), < * < n. 

Observe that the above expression coincides with Suetin-Rodrigues formula up to 
the binomial coefficients. 



Example 4: Orthogonal polynomials on the unit ball 

Classical orthogonal polynomials on the unit ball, B2 = {(x, y) : x 2 + y 2 < 1}, 
are associated with the weight function (symmetry factor) 

u J (x,y) = (l-x 2 ~y 2 r^ 2 , /i>-l/2. 
In this case, the matrices <E> and are given by 

r 1 xy \ ^_({2 i i + 2)x 



xy y 2 ~ir \{2n + 2)y 
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and the matrix Pearson-type equation for the symmetry factor to is given by 



2 



1 xy \ fuj x \ _ f(2[i-l)x 



xy y 2 - lj \uj y J \{2fi-i) y/ 

We can check that condition ((21(1 is not true, and matrix Rodrigues formula l(22|) 
does not provides a WOPS. However, we can transform the above matrix Pearson- 
type equation and obtain 

Equation 1)25(1 is a Pearson-type equation for the symmetry factor u/, it is diagonal, 
and it satisfies condition l(24|l . Then, matrix Rodrigues formula provides a WOPS 
defined as follows 

Qn = {Qn : 0, Qn— 1,1) Qo,n) , 

where 

(26) Q n _ i ^y) = ±( n ) d ^ d U(l- x 2-yi ru ), Q< l < n , 



1 1 

which coincides with the Rodrigues formula for the orthogonal polynomials on the 
unit ball as described in [H I1CJI fT5] . 

Example 5: Orthogonal polynomials on the simplex (Appell polynomi- 
als) 

The weight function associated with the classical orthogonal polynomials on the 
simplex, T — {{x, y) : x, y > 0, 1 — x — y > 0}, is defined by 

u{x,y) = x a y /3 (l - x - y) 7 , a,/3,7>-l. 

Here, the matrices $ and 'J are given by 

<T) ^/ar(x-l) xy \ $ _ ({a + (3 + 7 + 3)as - (a + 1) 



xy V(y-l)j' \(a + /3 + j + 3)y - {(3+ I) / 

Observe that condition (|21|l holds, and therefore, the matrix Rodrigues formula for 
triangle polynomials provides a WOPS. Since the obtained polynomials are monic, 
they coincide with the so-called second kind Appell polynomials (see |16|l 

Y. J. Kim et al (see JUJ) proved that to satisfies also a diagonal matrix Pearson- 
type equation 

'x (1 - x - y) \ ( lo x \ _ (a(l - x - y) - 7 x\ 

y(l-x-y)J \ujyj \f3(l - x - y) - jy) U ' 

In this way, condition l|21() holds, and we are again in the situation described in 
Example 1. Then, we get 

Qn = {Qn,0, Qn-1,1) ' ■ ■ ' Qo,n) j 

defined from the matrix Rodrigues formula 

(27) Q n _,, i (^y) = ^Qar i ^( a; n - i y i (l-x-yra;), < 1 < n, 
is a WOPS relative to u. 

Rodrigues formula for the circle and the triangle polynomials, l(26|l and 1(27(1 
respectively, coincide with the classical expressions for these polynomials which 
appear in the literature (see, for instance, [Tl HI UHl ITTI IT?)] ). 
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Example 6: The most intriguing case (sic, L. L. Littlejohn, |15| ) 
Krall and Sheffer showed that the differential equation 

L\p]= 3yp xx + 2p xy - xp x - yp y = -np, 

has an OPS as solutions. In this case, the matrices <E> and iff are given by 

.-(* i), 

A symmetry factor for the partial differential equation is ui(x,y) = exp(y 3 — xy) 
f|15|L Observe that condition (|21() is satisfied, and then we can obtain a matrix 
Rodrigues formula for these polynomials. In particular, we get 



j| I I IV..,, I •) 'J..-,, IV,.- c .. •-!...,- ■ (;„ 



1, 

6y,2a;y - 2,y 2 ) 
(-x 3 + I8xy - 12, -3x 2 y + 18y 2 + 6x, -3xy 2 + 6y, -y 3 ) . 



We must remark that Rodrigues-type formula for these polynomials can not be 
obtained using the Suetin tools (see JOj)- Moreover, the symmetry factor oj is not 
a weight function. These polynomials has attracted considerable attention fpiTTll 
EI El EH) since this is the simplest case of classical orthogonal polynomials in two 
variables orthogonal with respect to a non positive definite moment functional. 



7. Proof of the main result 

In this section, we will prove the matrix Rodrigues formula, introduced in The- 
orem m several steps that we will organize in a series of lemmas. 

Lemma 7.1. Let A — (dij)} j =0 be a 2 x 2 polynomial matrix. Assume that, there 
exist polynomial matrices *ff k = (4'i j)l j=o> such that 

(28) {a kfi A) x + {a k . 1 A) y =A$f k , k = 0, 1. 
Then, for n > 1, we have 

(29) (a M i W ) x + (a M 4 W ) y =i W ^, k = 0, 1, 

where ^ = (V , "j' C )^j=o ■M„+i('P), k = 0,1, are three-diagonal matrices with 

ipj-i,j = {n + l- j) ipo A , 1 < j < n, 

1>?f = {n-j)^lo+j^i,i-(n-l)[(akfi) x + (a k ,i) y ], < j < n, 
VftiJ = (7 + l)^,o. 0<j<»-1. 
Moreover, deg < deg A — 1, k = 0,1, n > 1. 

Proof. The Lemma follows by induction on n. For n = 1, the result holds using 
AW = A, and = Hf k , k = 0, 1. 

Let < i < n — 1, and < j < n, and assume that the result holds for n — 1. 
First, we compute the left hand side of (|29Jl using the induction hypothesis, and 
Recurrence I of Lemma 15.11 
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( a °.°4"0 x + ( a °' ia £ 



(a>o,o 



{n-l} {n-l} 
O0,0 O-i.j + a 0,l a i,j-l 



a>o,i 



{n-l} . {«-!} 
a 0,0 a ij + a 0,l a i,j-l 



0-0.0 

a>o,i 



(oo,o «ij + («o,i + [(ao,o)x ao,o + (ao,o)y ao,i] <4™ 1} 

[(ao,o a<ijZ{ 



An-l} 



)x + (&0,1 a i,j-l )y + [{ a 0,l)x ao,o + ( a 0,l)y a 0,l] 

a-o.o ^2 a l"~ 1} ^"j" 1 ' + { a °<° [^o,o + ( a o,o)x + (ao,i)»] + a o,i V'i.o} a l,j~ 1} 



3+1 



00,1 ^2 a M 1} + { a ofiiPo.i + ao,i [V 1 ?,! + (ao,o)s + ( a o,i) y ] } a|"_J } 

1=3-2 

{n-l} f ,n , /,n-l,0l . {n-l} /.n-1,0 

ao,o alj./ [^0,1 + Wj-x j \ + 00,1 V>j-2j-i 

^0,0 + ^jj 1 ' + ( a O,o)x + (OQ,l)i/ 

^1,1 + ipVltf-i + (ao,o)x + (ao,x)i/ 



{n-l} 
a 0,0 Ojj 



1 {n-l} 

+ ao,iai,j_i 

{n-l} , n-1,0 . {n-l} \ ,,0 , / n-1,0 

+ O0,0 Oij+l Vj+lJ + 00,1 Ojj ^1,0 + Vjj-l 

Now, replace the recurrence relations for ?/>™ _1 '°, and Recurrence I for the elements 
a} n ) to obtain 



(«o,o«g } ) x 



a i a 



{"} 



{n} , n,0 {n} , n,0 {«} , „,0 

a i,j-irj-ld^ a id +a ij + l^i+lj- 



The rest of the cases follows in a similar way. □ 

Lemma 7.2. Assume that the matrix polynomial $ satisfies \21)) . Let A n £ 

■ / ^(n+i)x(m+i)(^') ^ e arl arbitrary polynomial matrix. Then, 

div {n} (<!> {n} A„ u) = div {n - 1} (${"- x > A„_! w ), 
where A„_i G A / t nx (m+i)(^ : ') * s a polynomial matrix satisfying 

deg A„_i < dcgv4„ + 1. 

Proof. From the definition of div^ in J7J, we can easily deduce 

divW(${">A n w)=div^- 1 > [4(Ci$ W A l u) + fll ( (iii- 1 $ W A l w)' . 
Now, if we apply recurrence formulas in Lemma 15. II we get 
9 X (L°_! ^ A n w) + d y (L 1 n _ 1 A„ £j) = 
= "(a ${"- 1 >) ;c + (6 L°_ X A, w 

+ [(&<f^- 1 >) ;c + (c& n -V)y\ L l n _ x A n u 

[(od +bL 1 n _ 1 )(A n ) x + (6L°_! +cii_ 1 )(U] " 
+${«-!} ^(au^ + + £„_x An(6wi + cw„)] . 
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Condition J3TJ, Lemma \7. II and equation Ijl5(l give 

div { " } (${"} A n lo) = div { "~ 1} (pi*- 1 ) A„_ x w), 

where 

+(ai°_ x +bL 1 n _ 1 )(A n ) x + + C Li_ 1 )(4„) J . 
Observe that the degree condition can be deduced from the explicit expression for 

An-!. □ 

Now, taking A n — I n +i, and applying induction on n, we get 
Lemma 7.3. In the hypothesis of Lemma \7.ci\ for any n > 0, 

Qt n : = ^div { " } (${™>w), 

is a 1 x (n + 1) polynomial matrix of degree less than or equal to n. 

Using the same technique as Lemma l7.3l we can write in terms of the moment 
functional u. 

Lemma 7.4. Let u be a classical moment functional, and let |ZJ) be the matrix 
partial differential equation associated with u. Assume that condition $21]il holds. 
Then, for n>0, Q„ defined in Lemma \ 7. S\ satisfies 
(i) Q^u = div M (${">?/), n>0, 

(it) (u, QnX*,) =0, < m < n - 1, where X l m = (a;" 1 ^ y*)Ji „, 

fmj {Q^}« is a 1 X (n + 1) vector of polynomials of degree less than or equal to n. 
In fact, if we denote 

Qn — (Qn.Oj Qn-l,li ■ • ■ , Qo,n), 

then Qn-i.i has exact degree n or Q n -i,i = 0, for < i < n. 

Proof, (i) follows from a similar reasoning as used in Lemma 17.31 
In order to prove (ii), we compute 

(u,OA) = (Q n u,X t J = (div^(^u),X t J = 

= {-l) n (^ n K, V {n} X* m ) = 0, 0<m<n-l. 

As a consequence, (u,Q n p(x,y)) = 0, for any p(x,y) € V n -\- 

From the above property, we will prove (iii). Let n > 1, and suppose that there 
exists < i < n, such that deg Q n -i,i < n. Then, 

(u,Q n Q n -i,i) = 0, 

and therefore {u,Q' 2 n _ i 4 ) = 0. Since u is a quasi definite moment functional, we 
deduce Q n -i,i = 0. 

□ 



To end the Section, we will prove Corollaries 16.21 and 16.31 
Proof of Corollary 16.21 
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Proof. If $ is a diagonal matrix, 

Qn — (Qn,0, Qn-1,1, ■ ■ ■ , Qo,n), 

where 



Q n - k ,k(x, y)u = (J d n x - k d k y (a n ' k c k u), < k < n. 

From the previous Lemma, Q n -k,k has exact degree n or Q n -k,k = 0, for < 
k < n. Assume that Q n -k,k = 0, then d k (a n ~ k c k u) — 0, and therefore a = 
or c = which contradicts (|10fl . Moreover, 

(u, Q n - k , k (x, y) x n ~ l y l ) - (-1)" (fj {a n ~ k c k u, d^ k d k (x n ~ l y 1 )} = 0, 

if k 7^ I. In the case k = I, we have 

(u,Q n - k , k (x,y) x n - k y k )^0, 

since Q n -k,k has exact degree n. 

Finally, we are going to show that the polynomials {Qn-fc,fc,0 < k < n} are 
independent modulus V n -i- Let Ao, Ai, . . . , A„ be constants such that 

n 

q( x > v) = X] Afc Qn-kA x > y)> 

fe=0 

is a polynomial of degree less than or equal to n — 1. Since 

(u,q(x,y)p(x,y)) = 0, 

for any p € V n -\, we have (u, q(x, y) 2 ) = 0, and therefore q(x, y) = 0. 
Then, we get 

n 

= (u, q(x, y) x n ~ % y l ) = ^ A fc (Q n ^ k (x, y) u, x n ' 1 y l ) 

k=0 

= X i (Q n -i i i(x,y)u,x n ~ l y l ), 
so A, = 0, < i < n. □ 



Proof of Corollary EOl 

Proof. Observe that 

= u , V^ l >X^) = (-1)" n\ (u, 

and then, the result follows. □ 
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